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Abstract 

We propose and study an iterative substructuring method for an h-p 
Nitsche-type discretization, following the original approach introduced in [18] 
for conforming methods. We prove quasi-optimality with respect to the mesh 
size and the polynomial degree for the proposed preconditioner. Numerical 
experiments asses the performance of the preconditioner and verify the theory. 



1 Introduction 

Discontinuous Galerkin (DG) Interior Penalty (IP) methods were introduced in the 
late 70's for approximating elliptic problems. They were arising as a natural evo- 
lution or extension of Nitsche's method [15] . and were based on the observation 
that inter-element continuity could be attained by penalization; in the same spirit 
Dirichlet boundary conditions are weakly imposed for Nitsche's method [I5J. The 
use, study and application of DG IP methods was abandoned for a while, probably 
due to the fact that they were never proven to be more advantageous or efficient 
than their conforming relatives. The lack of optimal and efficient solvers for the 
resulting linear systems, at that time, surely was also contributing to that situation. 
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However, over the last 10-15 years, there has been a considerable interest in 
the development and understanding of DG methods for elliptic problems (see, for 
instance, [7] and the references therein), partly due to the simplicity with which the 
DG methods handle non-matching grids and allow for the design of hp-refinement 
strategies. The IP and Nitsche approaches have also found some new applications; 
in the design of new conforming and non- conforming methods [HI El EH1 123 EH EB] 
and as a way to deal with non- matching grids for domain decomposition [HI 153] . 

This has also motivated the interest in developing efficient solvers for DG meth- 
ods. In particular, additive Schwarz methods are considered and analyzed in [351 
[201 El SI El U US 120] . Multigrid methods are studied in [37J M\ M\ ■ Two-level meth- 
ods and multi-level methods are presented in [23 [22] and other subspace correction 
methods are considered in [121 [TU1 [TT] . 

Still the development of preconditioners for DG methods based on Domain De- 
composition (DD) has been mostly limited to classical Schwarz methods. Research 
towards more sophisticated non-overlapping DD preconditioners, such as the BPS 
(Bramble Pasciak Schatz), Neuman-Neuman, BDDC, FETI or FETI-DP is now at 
its inception. Non-overlapping DD methods typically refer to methods defined on a 
decomposition of a domain made up of a collection of mutually disjoint subdomains, 
generally called substructures. These family of methods are obviously well suited 
for parallel computations and furthermore, for several problems (like problems with 
jump coefficients) they offer some advantages over their relative overlapping meth- 
ods, and have already proved their usefulness. Roughly speaking, these methods 
are algorithms for preconditioning the Schur complement with respect to the un- 
knowns on the skeleton of the subdomain partition. They are generally referred 
sub structuring preconditioners. 

While the theory for the confoming case is now well established and understood 
for many problems [5T], the discontinuous nature of the finite element spaces at the 
interface of the substructures (in the case of Nitsche-type methods) or even within 
the skeleton of the domain partition, poses extra difficulties in the analysis which 
preclude from having a straight extension of such theory. Mainly, unlike in the 
conforming case, the coupling of the unknowns along the interface does not allow 
for splitting the global bilinear form as a sum of local bilinear forms associated 
to the substructures (see for instance [35] and [31 Proposition 3.2]). Moreover the 
discontinuity of the finite element space makes the use of standard iJ^-norms in 
the analysis of the discrete harmonic functions difficult. 

For Nitsche-type methods, a new definition of discrete harmonic function has 
been introduced in [29] together with some tools (similar to those used in the analysis 
of mortar preconditioners) that allow them to adapt and extend the general theory 
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[5Tj for substructuring preconditioners in two dimensions. More precisely, in [221 E01 
[3T| the authors introduced and analyzed Balancing Domain with Constrains BDDC, 
Neuman-Neuman and FETI-DP domain decomposition preconditioners for a first 
order Nitsche type discretization of an elliptic problem with jumping coefficients. For 
the discretization, a symmetric IP DG scheme is used (only) on the skeleton of the 
subdomain partition, while piecewise linear conforming approximation is used in the 
interior of the subdomains. In these works, the authors prove quasi-optimality with 
respect to the mesh-size and optimality with respect to the jump in the coefficient. 
They also address the case of non-conforming meshes. 

More recently, several BDDC preconditioners have been introduced and analyzed 
for some full DG discretizations [3H1 [221 E], following a different path. In [IS] 
the authors consider the p-version of the preconditioner for an Hibridized IP DG 
method [2U [M], for which the unknown is defined directly on the skeleton of the 
partition. They prove cubic logarithmic growth on the polynomial degree but also 
show numerically that the results are not sharp. The IP DG and the IP-spectral DG 
methods for an elliptic problem with jumping coefficient are considered in [U] and 
[23] , respectively. In both works the approach for the analysis differs considerably 
from the one taken in [2§l I3"U1 I3"T] and relies on suitable space decomposition of 
the global DG space; using either nonconforming or conforming subspaces. This 
allow the authors to adapt the classical theories for analyzing the resulting BDDC 
precondit ioners . 

In this work, we focus on the original substructuring approach introduced in [TSJ 
for conforming discretization of two dimensional problems and in [TjJJ [32] for three 
dimensions (see also [521 EI] f° r a detailed description). In the framework of non 
conforming domain decomposition methods, this kind of preconditioner has been 
applied to the mortar method [H [161 HS1 SB] and to the three fields domain decom- 
position method [U] , always considering the h- version of the methods. For spectral 
discretizations and the p version of conforming approximations the preconditioner 
has been studied in [HI H2j . For h-p conforming discretizations of two dimensional 
problems the BPS preconditioner is studied in [2]. To the best of our knowledge, 
this preconditioner has not been considered for Nitsche or DG methods before. 

Here, we propose a BPS (Bramble-Pasciak-Schatz) preconditioner for an h-p 
Nitsche type discretization of elliptic problems. In our analysis, we use some of the 
tools introduced in [211 [30], such as their definition of the discrete harmonic lifting 
that allows for defining the discrete Steklov-Poincare operator associated to the 
Nitsche-type method. However, our construction of the preconditioners is guided by 
the definition of a suitable norm on the skeleton of the subdomain partition, that 
scales like an if 1 / 2 -norm and captures the energy of the DG functions on the skeleton. 
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This allow us to provide a much simpler analysis, proving quasi-optimality with 
respect to the mesh size and the polynomial degree for the proposed preconditioners. 
Furthermore, we demonstrate that unlike what happens in the conforming case, to 
ensure quasi-optimality of the preconditioners a block diagonal structure that de- 
couples completely the edge and vertex degrees of freedom on the skeleton is not 
possible; this is due to the presence of the penalty term which is needed to deal with 
the discontinuity. We show however that the implementation of the preconditioner 
can be done efficiently and that it performs in agreement with the theory. 

The rest of the paper is organized as follows. The basic notation, functional 
setting and the description of the Nitsche-type method are given in next section; 
Section [2J Some technical tools required in the construction and analysis of the 
proposed preconditioners are revised in Section[3J The substructuring preconditioner 
is introduced and analyzed in Section HJ Its practical implementation together with 
some variants of the preconditioner are discussed in Section |5j The theory is verified 
through several numerical experiments presented in Section [6j The proofs of some 
technical lemmas used in our analysis are reported in the Appendix |A1 

2 Nitsche methods and Basic Notation 

In this section, we introduce the basic notation, the functional setting and the 
Nitsche discretization. 

To ease the presentation we restrict ourselves to the following model problem. 
Let Q C M 2 a bounded polygonal domain, let / G L 2 (Q) and let 



The above problem admits the following weak formulation: find u* e Hq(Q) such 
that: 




u* = f 
u* = 



in f2, 
on <9f2. 



a(u*,v) = f(v) 



for all v E H%(Q), 



(1) 



where 
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2.1 Partitions 



We now introduce the different partitions needed in our work. We denote by Th a 
subdomain partition of fi into N non-overlapping shape-regular triangular or quadri- 
lateral subdomains: 

N 

We set 

H £ = min ifg,- where Hgj = \dQe n dQj\ , (2) 

and we also assume that Hi ~ diam(fi^) for each £ = 1, . . . N. We finally define the 
granularity of Th by H = min^ Hi. We denote by T and T 9 respectively the interior 
and the boundary portions of the skeleton of the subdomain partition Th'- 

N 

T = |Jr £ , T e = dQ e ^dQ V£ = 1,...,N. 

i=i 

N 

r 9 = |Jrf, r 9 £ =:dQ e ndQ ve = i...,N. 

We also define the complete skeleton as £ = T U T 9 . The edges of the subdomain 
partition that form the skeleton will be denoted by E C T and we will refer to them 
as macro edges, if they do not allude to a particular subdomain or subdomain edges, 
when they do refer to a particular subdomain. 

For each fig, let {7^} be a family of /ine partitions of fi^ into elements (triangles 
or quadrilaterals) K with diameter Jik- All partitions 7^f are assumed to be shape- 
regular and we define a global partition Th of fi as 

Observe that by construction Th is a fine partition of fi which is compatible 
within each subdomain fi^ but which may be non matching across the skeleton T. 
Throughout the paper, we always assume that the following bounded local variation 
property holds: for any pair of neighboring elements K + G T h e+ and K~ G 7jf , 
£+^£-,h K+ ~h K -. 
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Note that the restriction of Th to the skeleton T induces a partition of each 
subdomain edge E C T. We define the set of element edges on the skeleton T and 
on the boundary of O as follows: 

£° h -. = {e = dK + n dK- n r, K + e if , K~ e Tf , £ + ^ r } , 
^: = { e = Mn90,Ke 7~4 , 

and we set £h = £%U£®. When referring to a particular subdomain, say O, for some 
£, the set of element edges are denoted by 

£f = {e E E° h : e C 50,}, = {e e £ 9 h : e C 90,}, ^ = ^ U S 3 / . 

2.2 Basic Functional setting 

For s > 1, we define the broken Sobolev space 

H s (7~h) = j0G L 2 (0) : <^e# s (0,) VO, eTffj-n^^)' 
whereas the trace space associated to H 1 (Th) is defined by 

For it = (u £ )f =1 in if 1 (7^) we will denote by U| E the unique element (ft = (cj/)f =l in 
$ such that 

We now recall the definition of some trace operators following [7], and introduce 
the different discrete spaces that will be used in the paper. 

Let e € £% be an edge on the interior skeleton shared by two elements K + and 
K~ with outward unit normal vectors n + and n~, respectively. For scalar and 
vector- valued functions ip E H 1 (Th) and r E [if 1 (Tff)] 2 , we define the average and 
the jump on e E as 

{r} = -(t+ + t-), [<p] = <p + n+ + <p-n- , on e E E° h 
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On a boundary element edge e G we set {r} = r and [<£>]= y^n, n denoting the 
outward unit normal vector to Q. 

To each element K G 7/jf, we associate a polynomial approximation order px > 1, 
and define the /zp-finite element space of piecewise polynomials as 

X l h = {ve C°(Q e ) such that v\ K e ¥ PK (K), Kerf}, 

where F PK (K) stands for the space of polynomials of degree at most px on K. 
We also assume that the polynomial approximation order satisfies a local bounded 
variation property: for any pair of elements K + and K~ sharing an edge e G 
Pk+ ^Pk-- 

Our global approximation space Xh is then defined as 

N 

X h = {v G L 2 (Q) : such that % G X e h } ~ J]*« , . 

We also define X° C -A^ as the subspace of functions of X h vanishing on the skeleton 
E, i.e., 

Xl = {v G X h : such that v\ x = 0}. 
The trace spaces associated to X[ and Xh are defined as follows 

= {7/ G H 1/2 (dn e ) : r/ = tu^ for some 10 G X^} W = 1, . . . , N 

N 

Notice that the functions in the above finite element spaces are conforming in the 
interior of each subdomain but are double-valued on T. Moreover, any function 
v G Xh can be represented as v = (v )f =1 with v G X g h . 

Next, for each subdomain Qe G Th and for each subdomain edge E C dQe, we 
define the discrete trace spaces 

= $i u , = {r/ G : 7/ = on aE}. 

Note that, since we are in two dimensions, the boundary of a subdomain edge E is 
the set of the two endpoints (or vertices) of E, that is if E = (a, b) then dE = {a, b}. 
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Finally, we introduce a suitable coarse space C $, that will be required for 
the definition of the subtructuring preconditioner: 

Zh = {v = (v e ) e $ ; F\E), V£ c dSl t , Vfi< e T H } . (3) 



2.3 Nitsche-type methods 

In this section, we introduce the Nitsche-type method we consider for approximat- 
ing the model problem fll]). Here and in the following, to avoid the proliferation of 
constants, we will use the notation x < y to represent the inequality x < Cy, with 
C > independent of the mesh size, of the polynomial approximation order, and of 
the size and number of subdomains. Writing x ~ y will signify that there exists a 
constant C > such that C~ l x < y < Cx. 

We introduce the local mesh size function h G L°°(E) defined as 

[h K if x G OK n 90, 

n 0) = < ,± (4) 

\ mm{h K+ , h K -} if x G &ftT + n dK~ n T, if G 7^ , ^ r , 

and the local polynomial degree function p G L°°(E): 

p x if x G OK n 90, 

(5) 

max{p^+,p^-} if x G n dK~ n T, K ± G 7; £± , ^ + ^ C. 

Remark 2.1. ^4 different definition for the local mesh size function h and the local 
polynomial degree function p involving harmonic averages is sometimes used for the 
definition of Nitsche or DG methods 129$ . We point out that such a definition yields 
to functions h and p which are of the same order as the ones given in and (TJP, 
and therefore result in an equivalent method. 

We now define the following Nitsche-type discretization [50j [H] to approximate 
problem ([T]): find u* h G X h such that 

A h (u* h , v h ) = f(v h ) for all v h G X h , (6) 

where, for all m,d6 X^, .A/^-, •) is defined as 

N 



A h (u, v) = ^2 { Vn i • I v 1 ds 



(7) 



Here, a > is the penalty parameter that needs to be chosen a > a f° r some 
a > 1 large enough to ensure the coercivity of Ah(-, •)■ 

On Xh, we introduce the following semi-norms: 

N 

MU = E W \<e h = E IIP h " 1/2 [^H^(e)> (8) 

together with the natural induced norm by Ah(-, ■)'■ 

\H 2 A=K,T B + ^K,e h Vvel,,. (9) 

Following [50] (see also 0) it is easy to see the bilinear form Ah(-, ■) is continuous 
and coercive (provided a > «o ) with respect the norm (j9j, i.e., 

Continuitiy : \A h (u,v)\ < \\u\\ A \\v\\ A \/u,veX h 
Coercivity : A h (v,v) > \\v\\ 2 A Wv e X h . 

From now on we will always assume that a > ckq. Notice that the continuity and 
coercivity constants depend only on the shape regularity of l~h- 



3 Some technical tools 

We now revise some technical tools that will be required in the construction and 
analysis of the proposed preconditioners. 

We recall the local inverse inequalities (cf. [3D], for example): for any rj e W K (K) 
it holds 

\v\H r (e) ^ Vk h K r \n\H'{e), e C dK 

for all s,r with < s < r < 1. Using the above inequality for s = and space 
interpolation, it is easy to deduce that for a subdomain edge E C <9f2^ and for all 
s, r, < s < r < 1, for all 77 G | B it holds that 

< max i (plO^^k'fe) <rf (s_r) /ir r |^k s (s), (10) 
\v\Hr( d n e )< max (Pk^ 1 )^^!^^) ^rf h s f r \r)\ s , dnv (11) 
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where hi and pi refer to the minimum (resp. the maximum) of the restriction to 
dQe of the local mesh size h (resp. the local polynomial degree function p), that is, 

ht = min h(x) and pg = max p(x). (12) 

xean e nr xedn e nr 

We write conventionally 

Hp 2 (H e 



h 



max . (13) 



The next two results generalize [HI Lemma 3.2, 3.4 and 3.5] and [T5j Lemma 3.2] 
to the hp- version. The detailed proofs are reported in the Appendix lAl 

Lemma 3.1. Let r\ = {rf)f =1 G $h and let \ = (x e )eLi £ £h be such that x £ ( a ) 
if [a] at all vertices a of Viz, for all Qi G Th- Then 

Lemma 3.2. Lei £ G $^ suc/i i/iai £(a) = a// vertices a of Qg. Let Cl G 

H 1 /' 2 {dVLe) be linear on each subdomain edge of dVL^. Then, it holds 



£ ii^, f(£) 4 + io g (^))Vai^ 

where hi and pi are defined in ( fl2|) and Hi is defined as in (j2J) . 



(dn e ) 



3.1 Norms on $^ 

We now introduce a suitable norm on $^ that will suggest how to properly construct 
the preconditioner. The natural norm that we can define for all rj = (rjije G $h is: 

inf \\u\\ A , (14) 
u G X h 

U\T, = V 

where the inf is taken over all u G Xh that coincide with 7] along S. We recall that 
since on T both u and r) are double valued, the identity r\ = u\ s is to be intended as 
n = u\ . Although 0141) is the natural trace norm induced on $ by the norm 0, 
working with it might be difficult. 
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For this reason, we introduce another norm which will be easier to deal with 
and which, as we will show below, is equivalent to (I14p . The structure of the 
preconditioner proposed in this paper will be driven by this norm. We define: 

IMIL* = M^/w) + a Yl IIP h ~ 1/2 MHi 2 (e)- (15) 

The next result shows that the norms (TH|) and (15. 2p are indeed equivalent: 
Lemma 3.3. The following norm equivalence holds: 

IMk ~ IMk* ^ IMk V77 G $ fe 

Proof. We first prove that IMk,* < IMk- Let r] = (77.^=1,. ..,jv G $h and let 
u = (u^)f=i jy such that iti£ = 77 arbitrary. Thanks to the trace inequality, we have 

11 2 < 1 1 2 

\ u £\H 1 (ri e ) > 

and so, summing over all the subdomains Qg G we have 

Adding now the term a ^ eg £ h llph -1 ^ 2 ! 7 ?!!^^) to both sides, and recalling the def- 
inition of the norms (I14p and (15. 2 p we get the thesis thanks to the arbitrariness 
of u. 

We now prove that |Mk < |Mk,*- Given r\ = (^)^=i ) ... ) £ G $/,, let w<? G 
be the standard discrete harmonic lifting of 77 , for which the bound lu^jfimg) < 
|r7^|^i/2(9Q f ) holds (see e.g. [IB]) and let u = (^)^=i,...,l Summing over all the sub- 
domains Vti and adding the term aY2 e e£ h Hp* 1 1 llL 2 (e) we S e ^ 

IMk < NU ~ IMk,*- 

□ 

4 Substructuring preconditioners 

In this section we present the construction and analysis of a substructing precondi- 
tioner for the Nitsche method (E])-([7j). 
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The first step in the construction is to split the set of degrees of freedom into 
interior degrees of freedom (corresponding to basis functions identically vanishing on 
the skeleton) and degrees of freedom associated to the skeleton T of the subdomain 
partition. Then,the idea of the "substructuring" approach (see [18]) consists in 
further distinguishing two types among the degrees of freedom associated to T : 
edge degrees of freedom and vertex degrees of freedom. Therefore, any function 
u G Xh can be split as the sum of three suitably defined components: u = u° + up = 



We first show how to eliminate (or condensate) the interior degrees of freedom 
and introduce the discrete Steklov-Poincare operator associated to ([7]), acting on 
functions living on the skeleton of the subdomain partition. We then propose a 
preconditioner of substructuring type for the discrete Steklov-Poincare operator and 
provide the convergence analysis. 

4.1 Discrete Steklov-Poincare operator 

Following [29||30], we now introduce a discrete harmonic lifting that allows for defin- 
ing the discrete Steklov-Poincare operator associated to ([7]). We also show that such 
a discrete Steklov-Poincare operator defines a norm that is equivalent to the one de- 
fined in (CT . 

Let X® C Xh be the subspace of functions vanishing on the skeleton of the 
decomposition. Given any discrete function w G Xh, we can split it as the sum of an 
interior function w° G X® and a suitable discrete lifting of its trace. More precisely, 
following USED], we split 



where, for rj G Rhijf) G X h denotes the unique element of X h satisfying 



u° + u E + u v . 



w = w° + R h (w\x), 




Rh(v)\*=V, A h (R h (v),v h ) = Wv h eX° h . 
The following proposition is easy to prove (see [291 130])- 
Proposition 4.1. For rj = (rj e ) G <ly, the following identity holds: 



(16) 
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and w° G X[ n if^f^) fremg £/&e solution of 

</q £ Jdn e 

The space can be split as direct sums of an interior and a trace component, 
that is 

X h =X° h ®R h {$ h ). 

Using the above splitting, the definition of Rh(-) and the definition of Ah(-, ■)> it is 
not difficult to verify that, 

A h (w, v) = A h {w°, v°) + A h (R h (w ls ), R h (v ]s )) 

= a(w°, v°) + s(w\ s , v\ s ), \/w,v&X h 

where the discrete Steklov-Poincare operator s : $^ x $^ 1 is defined as 

s(C, v) = A h (R h (0, Rh(v)) Vt,V£$h. (17) 
We have the following result: 
Lemma 4.2. Let Rh be the discrete harmonic lifting defined in ffT6l) . Then, 

\\Rh(v)\\A - \\vh h ,* V77G$/». 

Proof. If we show that H-R/it^lU — IMI* h j then the thesis follows thanks to the 
equivalence of the norms shown in Lemma 13.31 First, we prove that ||-R/,,(?7) ||.a < 
ll^ll*^; let rj G <&h, then from the definition of the inf, we get that 

3u G Xh ■ u\s = rj such that \\u\\a < 2||r/||$ h . 

Then, we can write Rh(jj) = u + v with v G X®, and ffl6|) reads 

A h (v, w) = -A h (u, w) Ww G X° . 

Setting w = v G X® in the above equation, leads to 

A h {v,v) = -A h (u,v) . 

Then, using the coercivity and continuity of Ah(-, •) in the || ■ ||^-norm we find 

\\ V \\ 2 A ~ <Ah{v,v) = \A h (u,v)\ < ||u|U||u|U • 
Hence, \\v\\a < IMU> an d so this bound together with the triangle inequality gives 

\\Rh(v)\\A < \\u\\a + IMU ^ \\u\\a ^ IMk ■ 
The other inequality < H-R/i^lU follows from the trace theorem. □ 
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From the above result, the following result for the discrete Steklov-Poincare 
operator follows easily. 

Corollary 4.3. For all £ G <&h, it holds 

s{t,t)*U\\l h ,.- 

Proof. Let £ G &h then from the definition of s(-, •), the continuity and coercivity of 
Ah(-, •) and applying Lemma I4T21 we have 

□ 



4.2 The preconditioner 

Following the approach introduced in [IS], we now present the construction of a 
preconditioner for the discrete Steklov-Poincare operator given by s(-,-). We split 
the space of skeleton functions as the sum of vertex and edge functions. We start 
by observing that C We then introduce the space of edge functions $f C $/,, 
defined by 

$f = {r] G = at all vertex A of ^ G 7^} 

and we immediately get 

$ h = £^ © $f . (18) 
The preconditioner s(-, •) that we consider is built by introducing bilinear forms 

s E : $f x $f — s> R s y : £ H x £^ — ► R 

acting respectively on edge and vertex functions, satisfying 

E wi^ v^G$f, (19) 

s V (v V ,V V )^ E I^IU(^) V^g£h, (20) 

and we define s : $^ x $^ — >• R as 

«(»?, = S E (v E , C E ) + S V (V V , n + q(r), 0, (21) 
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where 

4(77,77) = aJ2 llpk- 1/2 M||! 2(e) Vr? G $ fc . (22) 

Finally, we can state the main theorem of the paper. 

Theorem 4.4. Let s(-, •) and s(-, •) fre £/ie bilinear forms defined in (ITT)) and (l2Tj) . 

respectively. Then, we have: 

(l + log (Hp 2 /h))~ 2 s(r),r)) < 3(77,77) < §(77,77) Vr? G . 

The proof of Theorem 14.41 follows the analogous proofs given in [T51 [T5] for 
conforming finite element approximation. We give it here for completeness. 

Proof. We start proving that 5(77,77) < 5(77,77). Let 77 G then, r] = r] v + r] E with 
77^ G $f and rj v G By using Corollary 14.31 as well as the properties (|T9l) - (l20l) 
of the edge and vertex bilinear forms, and (|2"2"|) of q(-, •), we get 

s(r},V)<\\r}\\l h ,*= E l^ + ^lW^+«EllP h " 1/2 NHI^) 

~ E 1^1 i/2,sn,+ E l^liA^ + ?( 7 7» 7 7) 
< s 25 ^, 77 s ) + s v ( V v , V v ) + 5(77, 7/), 

and hence 

s(v,v) < 5(77,77) V77G$ h , . 

We next prove the lower bound. We shall show that 

S(V,V) £ (1 + log {Hp 2 /h)) 2 s(rj,rj) V77 G $ h . (23) 

For 77 G $/!, we have rj = t] v + r] E with 77 s G $f and 77 y G Then, from the 
definition of s(-, •) we have 

5(77, 77) = s E (r] E , t] E ) + s v (r] v , t] V ) + q(rj, 77) 
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Appling Lemma [3.21 with x = V E an d (l — V , we obtain 

I 1 2 



T.cT„ E^flO. O CT",, \ \ * / / 



that is 

^,^)<(i+io g (^^)) 2 E i^iU(^)- 

To bound s y (77 V , 77 ), we apply Lemma \3. II with x — V and if = 77, and we get 
s V (v V ,V V )< E l^lW^)<(l + log(^P 2 /^)) E MSrvw 

and hence 

W,^)+^(^,^)<(i+i 0g (if^)) 2 E m*v W 

Adding now the term a^ e&£h llph" 1 ^ 2 ! 7 ?]!! 2 ^) to both sides and recalling the defi- 
nition of q(-, •) we have: 

£(77,77) = s E (r) E ,r) E ) + s v (r] v ,r] v ) + q(rj,rj) 

< (1 + log (Hpyh)) 2 ( e i*/ W) +«E w^- 1/2 ivi\\i He) 

= {l + log{Hp 2 /h)) 2 \\r ] \\l i ^. 

Finally, using the equivalence norm given in Corollary 14 .3[ we reach (123|) and the 
proof of the Theorem is completed. 

□ 

As a direct consequence of Theorem 14.41 we obtain the following estimate for the 
condition number of the preconditioned Schur complement. 

Corollary 4.5. Let S and P be the matrix representation of the bilinear forms s(-, •) 
and §(■,■), respectively. Then, the condition number o/P _1 S ; k(P _1 S) ; satisfies 



V- l S)< (l + log(£Tp 2 //i)) 2 . (24) 
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Unfortunately, the splitting (TTHj) of $/j is not orthogonal with respect to the s(-, •)- 
inner product given in (|2~Tj) . and therefore the preconditioner based on s(-, ■) is not 
block diagonal, in contrast to what happens in the full conforming case. Furthermore 
the off-diagonal blocks in the preconditioner cannot be dropped without loosing 
the quasi-optimality. The reason is the presence of the q(-, •) bilinear form in the 
definition f[2~Tj) . and the fact that the two components in the splitting (I18p of $/j 
scale differently in the semi- norm that <?(•,•) defines. In fact, it is possible to show 
that, if for some constant k{K), it holds 

h V \\l h ,,<<h)\\ V \\l h ^ \/ r] = r ] v + r ] E e$ h , (25) 

then such n(h) must verify k{H) > H/h, which implies that, if we were to use a 
fully block diagonal preconditioner based on the splitting (fT8j) of $h an estimate 
of the form (|23|) would no be longer true. In order to show this, consider linear 
finite elements on quasi uniform meshes with meshsize h in all subdomains, and let 
V = {v e )e be the function identically vanishing in all subdomains but one, say ttk, 
and let r] k be equal to 1 in a single vertex of Qk and zero at all other nodes. With 
this definition, we have = |?7 fc | on dQk and = on S \ dQ k . Then, by a 

direct calculation, and recalling the definition of the semi-norm | • \*£ h in (JSJ), we 
easily see that 



i i o 1 i 1/ 1 9 h 

but \V \*,s h --r- 



or equivalent ly 



q(v v ,v v )-T^, qiviv)- 1 - (26) 

Therefore the energy of coarse interpolant t] v exceeds that of r) by a factor of Hk/hk- 
Hence, bounding t] v alone in the || ■ ||$ fe! *-norm would result in an estimate of the 
type (E5D 

q(v V ,V V )<\\v V \\l h ,,<<h)q(rj,rj), (27) 
which in view of fT26|) would imply 

Remark 4.1. We point out that the lack of the block-diagonal structure of the pre- 
conditioner associated to s(-,-) defined in WT\) , will not affect its computational 
efficiency, see Section® 
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5 Realizing the preconditioner 



We start by deriving the matrix form of the discrete Steklov-Poincare operator s(-, •) 
defined in (ITT)) . We choose a Lagrangian nodal basis for the discrete space Xh, and we 
take care of numbering interior degrees of freedom first (grouped subdomain-wise), 
then edge degrees of freedom (grouped edge by edge and in such a way that the 
degrees of freedom corresponding to the common edge of two adjacent subdomains 
are ordered consecutively) , and finally the degrees of freedom corresponding to the 
vertices of the subdomains. We let n i; n e and n v be the number of interior, edge 
and vertex degrees of freedom, respectively, and set n = n e + n v Problem ([6]) is then 
reduced to looking for a vector u G K ni+n with u = (u i; u e , u v ) solution to a linear 
system of the following form 

Aii A ie A lv \ /uA /FA 
A£ A ee A ev u e = F e . 
Al Al A vv / \u v J \F V J 

Here, u ± G M ni (resp. F ± G M ni ) represents the unknown (resp. the right 
hand side) component associated to interior nodes. Analogously, u e ,F e G M ne and 
u v , F v G M nv are associated to edge and vertex nodes, respectively. We recall that 
for each vertex we have one degree of freedom for each of the subdomains sharing 
it. For each macro edge E, we will have two sets of nodes (some of them possibly 
physically coinciding) corresponding to the degrees of freedom of Qf^ (E) and of 
K~(E). 

As usual, we start by eliminating the interior degrees of freedom, to obtain the 
Schur complement system 

with 

c _ f-^ee — A^ e A ±i A ie A ev — A^ e A i:L A lv \ _ — A^ e A 1:L FA 
" \A£ - Af v A^A ie A vv - A^A^AiJ ' g " [f v - A^A^F J " 

The Schur complement S represents the matrix form of the Steklov-Poincare opera- 
tor s(-, •). Remark that in practice we do not need to actually assemble S but only 
to be able to compute its action on vectors. 

In order to implement the preconditioner introduced in the previous section 
we need to represent algebraically the splitting of the trace space given by (I18p . 
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As defined in ([3]), we consider the space of functions that are linear on each 
subdomain edge, and introduce the matrix representation of the injection of 
into More precisely, we let E = {xj, i = 1, . . . , n e , n e + 1, . . . , n e + n v } be the 
set of edge and vertex degrees of freedom. For any vertex degree of freedom Xj, 
j = n e + 1, . . . , n e + n v , let <fj(-) be the piecewise polynomial that is linear on each 
subdomain edge and that satisfies 

V?j(x fc ) = 5 jk j, k = n e + 1, . . . , n e + n v . 

The matrix R T G M nXnv realizing the linear interpolation of vertex values is then 
defined as 

R T (i,j - n e + 1) = <pj(xi), z = l,...,n, j = n e + 1, . . . , n e + n v . 
Next, we define a square matrix R T G M nXn as 




I e G M neXne being the identity matrix. Let now S be the matrix obtained after 
applying the change of basis corresponding to switching from the standard nodal 
basis to the basis related to the splitting (fT8|) . that is 

S = RSR T = ^ | ve j . (28) 

Our problem is then reduced to the solution of a transformed Schur complement 

system 

Su = g, (29) 

where u = R~ T u and g = Rg. 



The preconditioner P. The preconditioner P that we propose is obtained as matrix 
counterpart of f l2~l]) . In the literature it is possible to find different ways to build 
bilinear forms s E (-,-), s v (-, •) that satisfy ( TT9"j) and fl2TJ|) . respectively. The choice 
that we make here for defining s E (-, •) is the one proposed in [IB] and it is based on 

1 /2 

an equivalence result for the H q norm. We revise now its construction. Let /o(-) 
denote the discrete operator defined on $°(-E) associated to the finite-dimensional 
approximation of —d 2 /ds 2 on E. It is defined by: 

(Io<P,<I>)e = (</,</>% V0G$°(S), (30) 
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where the prime superscript refers, as usual, to the derivative d/ds with respect to 
the arc length s on E. Notice that, since Iq(-) is symmetric and positive definite, its 
square root can be defined. Furthermore, it can be shown that 

n n /,i/2 a/2 

see [18]. Then, we define 

W,£ B )=£ E^V,^ Vv E ,ee$°e(E). (31) 

n e eT H Ecdn e 

For t] E e <&1{E) we denote by r\ E its vector representation. Then, it can be verified 
that, for each subdomain edge E C dQe, we have (see [T7] pag. 1110 and [28] ) 

where K E = M^ /2 (M^ 1/2 R B M^ 1/2 ) 1 / 2 M^ /2 , and where M E and R E are the mass 
and stiffness matrices associated to the discretization of the operator — d 2 /ds 2 (in 
with homogeneous Dirichlet boundary conditions at the extrema a and b of 

E. 

Observe, that for each macro edge E shared by the subdomains Qe+ and f^-, 
Kg is a two by two block diagonal matrix of the form 




where K. E are the contributions from the subdomains f^± sharing the macro-edge 
E. As far as the vertex bilinear form s v (-, •) is concerned, we choose: 




where %{■) denotes the standard discrete harmonic lifting [T8|l52]. We observe that if 
the Qgs are rectangles, for t] e we have that Hfj] is the Q 1 (fi^) polynomial that 
coincides with r] at the four vertices of Qg. Computing s v (r) v ,£ v ) for r] , !; v G 
is therefore easy, since it is reduced to compute the local (associated to Vt^) stiffness 
matrix for Q^f^) polynomials. 
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Remark 5.1. A similar construction also holds for quadrilaterals which are affine 
images of the unit square, and for triangular domains. In fact, if Qg is a triangle 
then for r\ G £# we have that T-L e h r] e is the P 1 ^) function coinciding with rf at the 
three vertices of fig. If fie is the affine image of the unit square, we work by using 
the harmonic lifting on the reference element. 

The preconditioner P can then be written as: 

















\ 











































K E M 







V 











p 

A vv 


/ 



where for each macro edge Ei 



Q, 



(33) 



K, 



: k ^) 1/2 



In (133]) P vv is defined as the matrix counterpart of (132]) whereas Q = RQR T and 



Q 



o 
o 



o 



T 

E 2 V 









QEiV \ 

Qe 2 v 

QvV) J 



is the matrix counterpart of ff22l . Remark that, due to the structure of the off 
diagonal blocks of Q, P is low-rank perturbation of an invertible block diagonal 
matrix. The action of P _1 can therefore be easily computed, see e.g. [2B] sec. 2. 7.4, 
p. 83. 



The preconditioner P*. For comparison we introduce a preconditioner P* with 
the same block structure of P but with the elements of the non-zero blocks coincid- 
ing with the corresponding elements of S. We expect this preconditioner to be the 
best that can be done within the block structure that we want our preconditioner 
to have. In order to do so, we replace the S ee component of S with the matrix 
obtained by dropping all couplings between the degrees of freedom corresponding to 
nodes belonging to different macro edges, and use the resulting matrix as precon- 
ditioner. More precisely, for any subdomain edge of the subdomain partition, 
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k = 1,. . . , M, let Jfc G ]g n exn e De (ii a g ona i matrix that extract only the edge 
degrees of freedom belonging to the macro edge Ek, i.e., 

-r / • .\ f 1 if i = j and Xj G E k . . 
Jfc(^j) = < n i,3 = l,- -,n e . 

I (J otherwise 

Then, we define 

m 

Pee = ^ ^ Jfc See Jfe 

This provides our preconditioner 



k=l 



P s 

1 ee u ev 



(34) 



Building this preconditioner implies the need of assembling at least part of the Schur 
complement; this is quite expensive and therefore this preconditioner is not feasible 
in practical applications. 

Remark 5.2. Note that we cannot drop the coupling between edge and vertex points, 
i.e. we cannot eliminate the off-diagonal blocks Qf^y, Q^.y- Indeed, as already 
pointed out at the end of Section \4-2\ with the splitting [W\l of &h it is not possible to 
design a block diagonal preconditioner without losing quasi- optimality. In Section® 
we will present some computations that show that the preconditioner 



D 



Pee 

o sL 



(35) 



is not optimal. 



6 Numerical results 

In this section we present some numerical experiments to validate the performance 
of the proposed preconditioners. 

We set Q = (0, l) 2 , and consider a sequence of subdomain partitions made of 



N = 4 squares, £=1,2,..., cf. Figure 1(a) for £ = 1, 2, 3, 4. For a given subdomain 



partition, £ = 1, 2, . . ., we have tested our preconditioners on a sequence of nested 
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structured and unstructured triangular grids made of n = 2 * 4 r , r = 
Notice that the corresponding coarse and fine mesh sizes given by H 

y-(r+l/2) _ _ 



+ 1, 



1,2, 



and h 



+ 1, 



respectively. In Figure 



la 



we have 

reported the initial structured grids, on subdomains partitions made by iV = 4 s 
2,3,4, are reported. 



squares, s = 

levels of unstructured grids on a subdomain partition made of N 



Figure 1(b) shows the first four refinement 

4 squares. 
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(a) Initial structured grids. 




(b) First three refinement levels of unstructured triangular grids on a subdomain partition made of 
N = 4 squares. 

Figure 1: Top: initial structured grids on subdomains partitions made by N = 4 e 
squares, t = 1,2,3,4. Bottom: first four refinement levels of unstructured grids on 
a subdomain partition made of N = 4 squares. 



Throughout the section, we have solved the (preconditioned) linear system of 
equations by the Preconditioned Conjugate Gradient (PCG) method with a rel- 
ative tolerance set equal to 10~ 9 . The condition number of the (preconditioned) 
Schur complement matrix has been estimated within the PCG iteration by exploit- 
ing the analogies between the Lanczos technique and the PCG method (see [361 
Sects. 9.3, 10.2], for more details). Finally, we choose the source term in problem 
(PP) as f(x,y) = 1, and set the penalty parameter a equal to 10. 

We first present some computations that show the behavior of the condition 
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Figure 2: Condition number estimate of the Schur complement matrix S versus 
1/h on different subdomains partitions made by N = 4 e squares, i = 1,2,3,4,5. 
Structured (left) and unstructured (right) triangular grids, piecewise linear elements 
(p = l). 



number of the Schur complement matrix S, cf. (jSJ). In Figure [2] (log-log scale) we 
report, for different subdomains partitions made by N = 4 e squares, £ = 1,2, 3, 4, 5, 
the condition number estimate of the Schur complement matrix S, k(S), as a func- 
tion of the mesh-size 1/h. We clearly observe that k(S) increases linearly as the 
mesh size h goes to zero. 

Next, we consider the preconditioned linear system of equations 

P _1 Su = P x g, 

and test the performance of the preconditioners P and P* (cf. ([55]) and (151]) . respec- 
tively). Throughout the section, the action of the preconditioner has been computed 
with a direct solver. 



In the first set of experiments, we consider piecewise linear elements (p = 1), and 
compute the condition number estimates when varying the number of subdomains 
and the mesh size. Table [T] shows the condition number estimates increasing the 
number of subdomains N and the number of elements n of the fine mesh. In Table [T] 
we also report (between parenthesis) the ratio between the condition number of the 
preconditioned system and (1 + \og(H/h)) 2 (between parenthesis). These results 
have been obtained on a sequence of structured triangular grids as the ones shown 



in Figure 1(a) Results reported in Table [TJ (top) refers to the performance of the 



24 



preconditioner P, whereas the analogous results obtained with the preconditioner 
P* are shown in Tabled] (bottom). We have repeated the same set of experiments 



Preconditioner P 



N I 


n — > 


n = 128 


n = 512 


n = 2048 


n = 8192 


n = 32768 


N = 


16 


3.11 (0.74) 


4.88 (0.65) 


7.50 (0.64) 


10.84 (0.64) 


14.79 (0.64) 


N = 


64 




3.30 (0.79) 


5.25 (0.70) 


8.00 (0.68) 


11.42 (0.67) 


N = 


256 






3.35 (0.81) 


5.36 (0.72) 


8.16 (0.70) 


N = 


1024 








3.37 (0.81) 


5.39 (0.72) 


Preconditioner P* 






n = 128 


n = 512 


n = 2048 


n = 8192 


n = 32768 


N = 


16 


2.26 (0.54) 


4.04 (0.54) 


7.01 (0.60) 


11.00 (0.65) 


15.83 (0.68) 


N = 


64 




2.42 (0.58) 


4.49 (0.60) 


7.85 (0.67) 


12.28 (0.72) 


N = 


256 






2.47 (0.59) 


4.60 (0.62) 


8.07 (0.69) 


N = 


1024 








2.48 (0.60) 


4.63 (0.62) 



Table 1: Preconditioner P (top) and P* (bottom). Condition number esti- 
mates and ratio between the condition number of the preconditioned system and 
(1 + \og(H I h)) 2 (between parenthesis). Structured triangular grids, piecewise linear 
elements (p = 1). 



on the sequence of unstructured triangular grids (cf. Figure 1(b)). The computed 
results are shown in Figure [2J As before, between parenthesis we report ratio be- 
tween the condition number of the preconditioned system and (1 + \og(H/h)) 2 . As 
expected, a logarithmic growth is clearly observed for both preconditioner P and P*. 



Next, always with p = 1, we present some computations that show that the 
preconditioner P^ defined as in f !35p . i.e., the block-diagonal version of the precon- 
ditioner P*, is not optimal (cf. Remark 15. 2p . More precisely, in Table [3] we report 
the condition number estimate of the preconditioned system when decreasing H as 
well as h. Table |3] also shows (between parenthesis) the ratio between k(PdS) and 
HhT 1 . We can clearly observe that on both structured and unstructured mesh con- 
figurations, the ratio between k^P^S) and Hhr 1 remains substantially constant as 
H and h vary, indicating that the preconditioner P^ is not optimal. 
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Preconditioner P 







n = 128 


n = 512 


n = 2048 


n = 8192 


n = 32768 


N = 


16 


2.87 (0.69) 


4.69 (0.63) 


7.35 (0.63) 


10.68 (0.63) 


14.62 (0.63) 


N = 


64 




3.05 (0.73) 


5.01 (0.67) 


7.75 (0.66) 


11.13 (0.66) 


N = 


256 






3.09 (0.74) 


5.08 (0.68) 


7.89 (0.67) 


N = 


1024 








3.11 (0.75) 


5.11 (0.68) 








Preconditioner 










n = 128 


n = 512 


n = 2048 


n = 8192 


n = 32768 


N = 


16 


1.84 (0.44) 


3.24 (0.43) 


5.51 (0.47) 


8.44 (0.50) 


12.00 (0.52) 


N = 


64 




2.01 (0.48) 


3.77 (0.50) 


6.35 (0.54) 


9.76 (0.58) 


N = 


256 






2.04 (0.49) 


3.90 (0.52) 


6.58 (0.56) 


N = 


1024 








2.05 (0.49) 


3.93 (0.53) 



Table 2: Preconditioner P (top) and P* (bottom). Condition number esti- 
mates and ratio between the condition number of the preconditioned system and 
(1 + log(H/h)) 2 (between parenthesis). Unstructured triangular grids, piecewise 
linear elements (p = 1). 



26 



Structured triangular grids 



n = 128 n = 512 n = 2048 n = 8192 n = 32768 



N = 


16 


11.51 (4.07) 


23.19 (4.10) 


47.40 (4.19) 


95.21 (4.21) 


190.69 (4.21) 


N = 


64 




11.58 (4.09) 


23.03 (4.07) 


47.16 (4.17) 


95.02 (4.20) 


N = 


256 






11.55 (4.08) 


22.96 (4.06) 


47.12 (4.16) 


N = 


1024 








11.44 (4.04) 


22.88 (4.04) 


Unstructured triangular grids 






n = 128 


n = 512 


n = 2048 


n = 8192 


n = 32768 


N = 


16 


9.45 (3.34) 


18.63 (3.29) 


39.13 (3.46) 


75.38 (3.33) 


148.93 (3.29) 


N = 


64 




8.93 (3.16) 


18.30 (3.24) 


38.88 (3.44) 


78.82 (3.48) 


N = 


256 






8.80 (3.11) 


17.85 (3.15) 


38.59 (3.41) 


N = 


1024 








8.75 (3.10) 


17.64 (3.12) 



Table 3: Preconditioner P^. Condition number estimates and ratio between /^P^S) 
and Hh" 1 (between parenthesis). Structured (top) and unstructured (bottom) tri- 
angular grids, piecewise linear elements (p = 1). 



Finally, we present some computations obtained with high-order elements. As 
before, we consider a subdomain partition made oi N = 4 e squares, £ = 1,2,..., 



(cf. Figure 1(a) for t = 1,2,3). In this set of experiments, the subdomain partition 
coincides with the fine grid, i.e., H = h , and on each element we consider the space 
of polynomials of degree p = 2,3, 4, 5, 6 in each coordinate direction. 

Table H] shows the condition number estimate of the non-preconditioned Schur 
complement matrix and the CG iteration counts. We have run the same set of ex- 
periments employing the preconditioners P and P*, and the results are reported in 
Table |5] . We clearly observe that, as predicted, for a fixed mesh configuration the 
condition number of the preconditioned system grows logarithmically as the polyno- 
mial approximation degree increases. Comparing these results with the analogous 
ones reported in Table IU it can be inferred that both the preconditioners P and P* 
are efficient in reducing the condition number of the Schur complement matrix. 
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N = n 


p = 2 


p = 3 


p = 4 


p = 5 


p = 6 


4 


5.1e+l ( 5) 


2.7e+2 ( 8) 


6.2e+2 (13) 


1.4e+3 (18) 


3.4e+3 (28) 


16 


3.2e+2 (22) 


8.4e+2 (42) 


2.0e+3 (69) 


4.6e+3 (101) 


l.le+4 (153) 


64 


1.2e+3 (90) 


3.2e+3 (150) 


7.6e+3 (231) 


1.8e+4 (312) 


4.3e+4 (446) 


256 


4.7e+3 (195) 


1.3e+4 (294) 


3.0e+4 (462) 


7.0e+4 (634) 


1.7e+5 (886) 



Table 4: Condition number estimates k(S) and CG iteration counts (between paren- 
thesis). Cartesian grids. 



Preconditioner P 



N = 


n 


p = 2 


p = 3 


p = 4 


p = 5 


p = 6 


N = 


1 


7.14 (1.25) 


9.04 (0.88) 


12.06 (0.85) 


14.15 (0.79) 


16.48 (0.78) 


N = 


16 


9.24 (1.62) 


9.93 (0.97) 


15.25 (1.07) 


15.99 (0.90) 


20.25 (0.96) 


N = 


64 


10.03 (1.76) 


10.14 (0.99) 


16.34 (1.15) 


16.57 (0.93) 


21.53 (1.02) 


N = 


256 


10.24 (1.80) 


10.19 (1.00) 


16.61 (1.17) 


16.71 (0.94) 


21.84 (1.04) 


Preconditioner P* 


N = 


n 


p = 2 


p = 3 


p = 4 


p = 5 


p = 6 


N = 


4 


1.88 (0.33) 


2.56 (0.25) 


3.75 (0.26) 


4.64 (0.26) 


5.70 (0.27) 


N = 


16 


4.60 (0.81) 


5.23 (0.51) 


8.71 (0.61) 


9.38 (0.53) 


12.25 (0.58) 


N = 


64 


6.18 (1.09) 


6.03 (0.59) 


10.35 (0.73) 


10.79 (0.61) 


14.33 (0.68) 


N = 


256 


6.55 (1.15) 


6.25 (0.61) 


10.83 (0.76) 


11.20 (0.63) 


14.94 (0.71) 



Table 5: Preconditioner P (top), P* (bottom). Condition number estimates and 
ratio between the condition number of the preconditioned system and (1 + log(p 2 )) 2 
(between parenthesis). Cartesian grids. 
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A Appendix 

In this section, we report the proofs of Lemma 13.11 and Lemma 13.21 



In the following, for E C subdomain edge we will make explicit use of the 
space Hq(E), < s < 1/2, which is defined as the subspace of those functions r\ 
of H S (E) such that the function fj G L 2 (dQg) defined as fj = t] on E and fj = on 
dfli \ E belongs to H s (dVLi). The space Hq(E) is endowed with the norms 

IMIffgCfO = ||^||H s (an«)- 

We recall that for s < 1/2 the spaces H S (E) and Hq(E) coincide as sets and have 
equivalent norms. However, the constant in the norm equivalence goes to infinity as 
s tends to 1/2. In particular on the reference segment E = (0, 1), for all </? G H S (E) 
and for all G R, the following bound can be shown (see |15j ) 

\v\m(E) ~ 1 ,1 z \\f - P\\h'<& + 



which, provided <p G H 1 ^ 2 (E), implies the bound 

1 1 

Prior to give the proofs of Lemmas 13.11 and 13.21 we start by observing that the 
following result, that corresponds to the hp- version of pJ3 Lemma 3.1], holds. 

Lemma A.l. Let E = (a,b) be a subdomain edge ofQg. Then, for all t] G $>e(E), 
the following bounds hold: 

(i) 

(V(a) -V(b)) 2 < (l + log (^)) MJrvw ( 37 ) 

(ii) If rj(x) = at some x G E it holds 

Nl!^)< (l + log (^)) M*v W (38) 

(iii) if r] £ §®(E), we have 
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Proof. We first show (ii). Notice that since E is an arbitrary subdomain edge, 
E C dQi, and so \E\ ~ H £ . We claim that for any <p G H 1 ^ 2+£ (E) the following 
inequality holds: 

I Ml £<»(£) < ^ l^ll^i?) + — |y|frVa+«(B)- ( 40 ) 

To show poj) one needs to trace the constants in the Sobolev imbedding between 
H 1 ^ 2+£ (E) and L°°(E). Let i£ =]0, 1[ be the reference unit segment. Then, for 
any if G H l ^ 2+£ (E), the continuity constant of the injection H^ 2+e (E) C L°°(E) 
depends on e as follows (see pJ3, Appendix], for details) 

ii - i|2 < || - ||2 | _| - |2 

A scaling argument using \E\ ~ Hi leads to (140|) . 

Let now i| e $/, and /3 G M an arbitrary constant. Using the inverse inequality 
(jTUI) . we have 

rr2 £ tt2e zr2e 4eh~2e 

1 In - e U 2 < £ Pi £ \n\ 2 

/»■ 2x (41) 

( H m\\ |2 

where in the last step we have taken e = 1/ \og(Hff> 2 jht) and used the fact that 
s V lc, g( s ) = e . Applying now inequality fHOj) to (f = 7] — (3 together with the above 
estimate (jUJ) yields: 



ll»? " PWUe) < Hi% ~ P\\h(B) + log (^) lCv W (42) 

Following [18], let /3 be the average over E of rj (or the L 2 - projection onto the space 
¥°(E)) of constants functions over E. Poincare-Friederichs inequality (or standard 
approximation results) give 

H t \\V ~ (3\\l*(e) <\v\m/*(E) (43) 

which yields 

\\V-P\\ 2 L-W< + (44) 



30 



The proof of (I38p is concluded by noticing that if r](x) = for some x G E then it 
follows 

|/3|<lh (45) 

which yields ( 138|) using triangular inequality 

The proof of (|37|) follows by applying the estimate (|38|) to the function 77 — 77(a), 
which by hypothesis vanishes at a £ £. 

To show (m^, we first notice that for rj G <&1(E), we can always construct an 
extension rj Q such that 

Vo = Vo OTI - E rj — on dft^ \ i? . 

Using now the inverse inequality (jlip . we obtain the following bounds 

IMIif^ a (£) ~ \Vo\m/z(dn e ) < P 2 i £h I £ \Vo\m/2-z(dn e ) ^ P/hi £ \Vo\ H ^-^ E y (46) 

where the second inequality follows from the boundedness from Hq(^(E) to ^^{dVli) 
of the extension by 0. 

To estimate now the Hl /2 - £ {E) seminorm of t] we observe that fl36l) rescales as 

\^h^(e) Z ^ (V /2 llv? - £IU^) + W - P\ H ^ {E) ) + 

Taking now if = r] and choosing /3 as its average on E, the first term on the right 
hand side above is bounded by means of Poincare-Friederichs inequality, and the 
second by means of estimate (j45]) . which holds since rj (a) = 0. Hence, we get 

.HfpfhJ e . HfvfhJ e ,, _.. 

\\^WhI' 2 (E) ~ ; \Vo\hV2(E) + /= \\Vo - P\\L°°(E)- 



Arguing as before and taking e = 1/ \og{Hep1/hi), and using bound ( 1441) we obtain 

II II < fl.Ll I I 

\Wo\\hI(?{E) ~ l^ 1 + l0g_ ^J I^UV2 (E) . 

Finally, since Y,Ecdn t I " \hV*(e) ~ I ' IjjV^oy b y squaring and taking the sum over 
E C 0ft*, we obtain ([39}. □ 



We are now able to prove Lemma 13.11 and Lemma 13.21 
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Proof of Lemma VJ. 11 A direct computation using the linearity of x shows that, if 
(X{ « bi are the vertices of the 2-th sub domain edge E % of Qg, we have 



i=l 

with N l E {= 3 or 4) denoting the number of subdomain edges of fig. Now, using f[3"TI) 
and assembling all the contributions we easily conclude that the thesis holds. 

□ 

Proof of Lemma \3.B. Let Co 3>£ be the unique element of $^ satisfying Co(°0 — 
for all vertices a of Vtg and (Co, r) H i/2 {dni) = (Cl, r) H i/2 {m ^ for all r G ¥ h with 
r(a) = for all vertices a of fig. It is not difficult to see that | ■ \H 1 / 2 (dQi) * s a norm 
on the subspace of functions in $^ vanishing at the vertices of fig and then, by 
standard arguments we get that Co is well defined and \Co\H 1 / 2 {an e ) ~ ICd/p/a^^). 
Now we can write: 

n% Ni N% 

E ^ {Eil S E »f + + E IKolli;,;^.,. < 47 > 

j=l i=l j=l 

with iV| number of subdomain edges of fig. The first sum on the right hand side of 
(1471) can be bound by using the previous lemma as 



Eii^c»f H , f(£;) <(i + io g (^i)) ! |« + aii, 
<( 1 + log (^!))Vaii 



/ 2 (£>n f ) 



where on one hand we used Poincare inequality to bound the if 1//2 norm of £ + Co 
(which vanishes at the vertices of fig) by the corresponding seminorm, while the last 
inequality follows by observing that, by the definition of Co, £ + Co £ $1 vanishes at 
the vertices of Qg and satisfies (Cl — Co,£ + Co^i/a^f) = 0- Hence, we have 

\C + C-L|ffi/2(ao^ = |£ + Co|^i/2( 9 q^ + |Cl _ Co|ffi/2(an^) - l£ + Co I ^1/2^) • 

Let us now bound the second sum on the right hand side of (T4T1) : we first observe 
that 

H^llfl^CE*) = I^°IhV2(£;|) + ^i(Co) + MCo), 
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having set 

with aj and 6, the two vertices of the subdomain edge E\. Now we can write 

^ KoIhV2(e|) < ICoIhi/z^) ~ 10^1/2(0^) < XX L ( a *) ~~ Cl(^)) 2 

i=l i=l 

^(i+iog(^))if + au (3n() , 

where the inequality iCd^i^^) < Z)i^(Ci(°i) ~ Chih)) 2 is proven in [18] by direct 
computation, and the last inequality follows by applying the bound of Lemma IA.1I - 
(ii) to the function (£ + — (£ + Cl)(^)- 

Let us now bound Ji. For notational simplicity let us identify Oj = and &j = i/. 
Adding and subtracting Cl^) + Cl(0) an d using the Cauchy-Schwarz inequality, we 
have 

JlKo) = /" «^dx< ric.w-aM + a(o)P dx+ ^iaw-a(o)P dx 

</o p| </o PI Jo \ x \ 

(48) 

Let us bound the first integral on the right hand side of (1481) . Setting Cl = Co — (l, 
we have 

H \c±(x) - c±m\ r h law - am 2 , , ic±(*) - 



dx = / : : dx + / : : dx. 



X Jo UM ./ /, | Jb | 

The first term can be bounded by 
/* ia(s)-Cx(0)l a /" N/otCJxMdrr 

7o N ~ 7o M ^K-J-Ihi^) ~ ICj-lffi/a(Bi 

while we bound the second term by 

A h dx^ia-awii^iog^^-j 



£ log -^ K 



2 



|2 



V he J J '^'^W 
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Next, we estimate the second integral on the right hand side of ( 1481) . By direct 
calculation and using the linearity of (l, we have 

f H \(l{x) -Cl(0)| 2 , <( f (n s sfA\\2<r-i ( H tPt \ \t , t\* 
J b] ~ \Cl{Bi) - (l{Ai)) < log J \Cl + £\ B i/2 {E iy 

Hence, we conclude that 



^(Co)< (i + iog(^)) 2 Ko-ai^) + iog(^) ia+ei^ w 

The term I2 can be bounded by the same argument. Collecting all the previous 
estimates the thesis follows. □ 
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